Long range statistical fluctuations of the crossed Josephson current 
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We investigate the crossed Josephson effect in a geometry consisting of a double ferromagnetic 
bridge between two superconductors, with tunnel interfaces. The crossed Josephson current vanishes 
on average because the Andreev reflected hole does not follow the same sequence of impurities as 
the incoming electron. We show that i) the root mean square of the crossed Josephson current 
distribution is proportional to the square root of the junction area; and ii) the coherent coupling 
mediated by fluctuations is "long range" since it decays over the ferromagnet phase coherence length 
l v , larger than the exchange length. We predict a crossed Josephson current due to fluctuations if 
the length of the ferromagnets is smaller than and larger than the exchange length 

PACS numbers: 74.50.+r, 74.78. Na,74.78.Fk 

I. INTRODUCTION 

Transport properties of hybrid structures consisting of a superconductor (S) multiply connected to several normal 
metal (N) or ferromagnetic (F) electrodes has focused an important interest recentljiS^. In usual Andreev reflection 
at a single NS interface, a spin-up electron coming from the N side is reflected as a hole in the spin-down band 
while a Cooper pair is transferred in the superconductor. Multiterminal structures allow "non local" processes, in 
which a spin-up electron in one electrode is Andreev reflected as a hole in the spin-down band in another electrode, 
corresponding to non local transmission in the electron-hole channel. Conversely non local transmission in the electron- 
electron channel corresponds to a process in which a spin-up electron from one electrode is transmitted as a spin-up 
electron in another electrode. Transport theory of three-terminal FSF junctions including non local transmission 
in the electron-electron and electron-hole channels has been discussed recently^SiSi^iiliiiiiiSiiii^iiSiifiiiLiS, in the 
tunnel linrft^ and for highly transparent interfaces^, on the basis of microscopic Green's functions^SiLSiSiiSi and 
in the framework of the scattering approachii*i2*ii The models were also extended to describe disorde»i2ii£*i&, non 
collinear ferromagnets^, and the nois o 17 ' 18 . On the experimental side, two experiments probing non local transport 
were carried out recentljii2i2£, in FSF and NSN three-terminal junctions. 

The question arises of whether the phase coherence of crossed Andreev reflection can be probed experimentally. 
We show here that this is possible with a double ferromagnetic bridge between two superconductors (see Fig. f^|. 
We have already shown that the crossed Josephson current vanishes on average in the diffusive limitS, because 
the Andreev reflected hole does not follow the same sequence of impurities as the incoming electron since they 
propagate in different electrodes. However, by evaluating the statistical fluctuations of the dc crossed Josephson 
current, we show here that the fluctuations of the Josephson current decay over the phase coherence length l v in the 
ferro magnet, larger than the decay length of the local average Josephson current set by the exchange length ^ (see 
Refs. !2ll22l23l24l25l2fl27l28l29l8(Tl . The fluctuations of the crossed supercurrent do not show 7r-shift oscillations and 
damping as a function of the ferromagnet length, because the spin-up and spin-down electrons of correlated pairs 
extracted from one superconductor do not see the same realization of disorder, so that the center of mass momentum of 
the spatially separated correlated pair averages to zero after propagation over a length comparable to the elastic mean 
free path. The root mean square of the crossed Josephson current is proportional to the square root of the junction 
area, because the number of diagrams involved in the supercurrent is equal to the junction area divided by the Fermi 
wave-length. The crossed Josephson supercurrent due to fluctuations can in principle be detected experimentally, 
provided the length of the ferromagnets is smaller than l v and larger than 

The article is organized as follows. Preliminaries regarding Green's functions are given in section^] The analytical 
results are presented in section ITTT1 for the average local supercurrent, and in section Hvl for the statistical fluctuations 
of the supercurrent. Concluding remarks are given in section [3 Some details on disorder averaging are provided in 
the Appendix. 



I 



* regis.melin@grenoble.cnrs.fr 

t U.P.R. 5001 du CNRS, Laboratoire conventionne avec l'Universite Joseph Fourier 



2 



Superconductor 



Ferromagnet 




sis 



Ferromagnet 




Superconductor 

FIG. 1: (Color online.) Schematic 3D representation of the Josephson junction considered in the article. 
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FIG. 2: Schematic 2D cut of the junction on Fig.0 We have represented some pairs of sites at the interfaces: (a^, ctk), (a[, a[) 
(6 m , /?m), and (b' n ,/3' n ). 



II. TECHNICAL PRELIMINARIES 



A. The models 



The superconductor is described by the BCS Hamiltonian 

„+ 



HbCS = X) _< ( C £<r C /V + 4 a C a ,a) + A X ( C ir C ii + C ",lCa,T 



(1) 



where i is the hopping amplitude, A is the superconducting gap, and a and f3 correspond to neighboring sites on 
a cubic lattice with a parameter ao. The lattice parameter ao is taken equal to the Fermi wave-length Xf- The 
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ferromagnetic electrodes are described by the Stoner model 



Tist 



(2) 



(a,/3),<7 



where /i ex is the exchange field. The exchange fields in the two ferromagnets are equal in the parallel alignment, and 
opposite in the antiparallel alignment. We incorporate also disorder scattering, described by the Hamiltonian 



'Hdh 



(3) 



where the impurities are located at the sites a n . The impurity scattering potentials V an are random variables. The 
site afc is on the ferromagnetic side of the interface, and the site on the superconducting side. 

The couplings between the ferromagnets and the superconductors are described by the tunnel Hamiltonian. The 
tunnel Hamiltonian at the interface (a, a) takes the form 



^a.a — ^ t [ ta k ,a k c aktCri F c a k ,cr,S ^a k ,a k c Qfc a ^g c a k ,o,F 



k.a 



(4) 



where the summation runs over all sites at the interface (see Fig.[2J, and where i afejQfc = t ak ,a k is the hopping amplitude 
connecting the sites ak and a^- 



B. Green's functions of a ferromagnet and a superconductor 



The starting point is the ballistic Green's function gi j (to) of the isolated ferromagnetic and superconducting elec- 
trodes in the Nambu representation. The ballistic Green's functions of a ferromagnet take the form 
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(5) 



(6) 



where g b and g b are the Green's functions of a spin-up electron and a hole in the spin-down band respectively, 
both having S z — 1/2, d a ^ is the distance between the sites a and b, ui the energy with respect to the chemical potential, 

T I T I 

Pf the density of states, k' F and k F the spin- up and spin-down Fermi wave- vectors, v F and v F the spin-up and spin- 
down Fermi velocities. The ballistic ferromagnet Green's functions given by Eqs. JSJ and J(|J) decay exponentially 



over the phase coherence length l^ all \ introduced phenomenologically through an imaginary part Hvf /Dp Ml ' to the 
energy lu. We note &f and vf the Fermi wave- vector and the Fermi velocity in the absence of spin polarization. We 
neglect in the following the energy dependence of the ferromagnet propagators in Eqs. JSJ and JfJJl since we suppose 
that the length R of the ferromagnets is small compared hv F /A and Hv F / A, both length scales being comparable to 
the ballistic BCS coherence length Hvf/ 'A. 

The Nambu Green's function of a ballistic isolated superconductor in the sector S z = 1/2 takes the form 
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where ps is the normal state density of states of the superconductor, d a fi the distance between the sites a and /3, and 



^bcs ( w ) = ^-f/v 7 A 2 — lu 2 the BCS coherence length at a finite energy. The information about propagation in the 
superconductor in the non local Josephson effect is contained in f a fi{to) = g^^Lu), where "1" and "2" refer to the 

electron and hole Nambu labels respectively. The statistical fluctuations of the Josephson current involve (f a fi{u)) 2 , 
where the overline is an average over disorder and over the different conduction channels. We have 32 ' 33 
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where ^bcs( w ) is the diffusive limit superconducting coherence length. 
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FIG. 3: (Color online.) Schematic representation of the lowest order diagrams for the "local" supercurrent (a), the "non local" 
supercurrent (b), and the statistical fluctuations of the "non local" supercurrent (c). We represent schematically the ladder 
diagrams for the diffuson. The rungs of the ladders correspond to a disorder scattering. 



C. Supercurrent 

The fully dressed Green's functions Gij (u>) are obtained from the Dyson equation G(u>) = g(w) + g(w) <8> £ ® G(u>), 
where ® corresponds to a summation over all the sites in the tunnel Hamiltonian The self-energy is provided by 
the couplings of the tunnel Hamiltonian, that, in the Nambu representation, take the form 



t a exp (iip/i) 

-t a exp(-i<p/4) 



(9) 



where tp is the phase difference between the two superconductors and t a is a real number. The phase does not 
evolve in time since we restrict here to the dc-Josephson effect. The second order diagrams on Fig. |3| contributing 
the supercurrent acquire a phase exp (Hip), giving rise to a supercurrent proportional to simp. The equilibrium 
supercurrent through electrode "a" is given by 
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(10) 



where the trace is a summation over the Nambu labels and the different conduction channels. The term (h cx — > —h cx ) 
corresponds to the contribution in the sector S z = —1/2. Eq. (|l(Jfl can be demonstrated from the expression of the 
current in terms of the Keldysh Green's function22i2£. The transparency of a single junction in the normal state 
is proportional to (i/e^) 2 , where is the Fermi energy. We suppose here that t <C €f, so that the supercurrent 
is expanded to order t A p 2 s p 2 F . The tunnel supercurrent coupling coherently the two ferromagnets is provided by the 
emission of a correlated pair of electrons from the left superconductor by Andreev reflection, followed by the absorption 
of the correlated pair by an Andreev reflection at the right superconductor. These two Andreev reflections can be 
"local" or "non local" in the sense that the incoming electron and outgoing hole can propagate in identical or in 
different electrodes. 

The local term 1^°^ in the supercurrent involves a diagram with propagation in a single ferromagnet (see Fig.|31-(a)), 
such that the incoming electron and the Andreev reflected hole are scattered by the same sequence of impurities: 
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(ii) 



where the overline is a disorder averaging, and where the sites a and b belongs to the left and right interfaces 
respectively (see Fig. To obtain Eq. (|ll|l . we start from Eq. (| 1 - use several times the Dyson equation, and 

replace the Green's functions of the isolated ferromagnets and superconductors by Eqs. (0), © and 10. The integral 
over energy is then calculated by contour integration. The poles of the products of the anomalous Green's functions 
are at uj — A, so that the ferromagnet Green's functions are evaluated at ui — A in Eq. 111|) . 

The spectral "non local" supercurrent involves a diagram with propagation in both ferromagnets (see Fig. 131(b) ): 



T-(nonloc) / \ 
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T,1,A/ \ T.2.A/ \ , T,2,A/ \ T,1.j4./ si 1 



(12) 
(13) 



where the sites a, 6, a' and b' belong to different interfaces (see Fig.|3-(b)). The variance of the non local supercurrent 
is obtained by integrating the square of the spectral supercurrent given by Eq. i|13f) over energy and averaging over 
disorder: 



(4 nonloc) ) = / ^(4 nonloc) H)' 



(14) 



III. EVALUATION OF THE AVERAGE LOCAL SUPERCURRENT 



The average over disorder of the product g1'1' A 9b'a' A m ^q- ^Hfl . is evaluated in the ladder approximation in the 
Appendix. We find^l 



1,1, A 2,2, A 



.'/,../, .'//,.., = ~ , 2 , T exp (iK h d a , b ) exp (-d a>b /£ h ), 



(15) 



where d a ,b is the distance between the sites "a" and "b" . The decay length of the supercurrent of a single SFS junction 
in the diffusive limit is given by^i 
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(16) 



where the exchange field enters through Afc = k F — kp, equal to the difference between the spin-up and spin-down 
Fermi wave vectors, with Afc = 2h/vp. The wave vector of the supercurrent oscillations£i*22i2&2ii25i2£i2Ii2£i22i2i is 

given by 



(ball) 



+ (Afc) 2 - 



,(ball) ' 



(17) 



where Id is the elastic mean free path. Due to the exchange field, the decay length £/, given by Eq. (jl6|) is smaller 
than the phase coherence length 



// /(ball) 
Ldi'V 



(18) 



Eqs. (|16fl and l|17|) follow from the identity 



— +iK h = 



(ball) 



iAk), 



(19) 



demonstrated in the Appendix. 

After summing over all pairs of sites (a, b) at the interface, as detailed in the Appendix, we obtain 



4 Z Z i 
rO°c) o^ E at Alt l 2 l+ I 2 - PsPF [ 

h k F l d a ^(1/^)2 + K 2 



(20) 



exp (-R/£ h ) cos (iT/j-R 



with tan6> = Kh/£,h- The local supercurrent is then proportional to the number of channels N c h at a single interface, 
7V c h being itself proportional to the area of a single junction. 
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IV. STATISTICAL FLUCTUATIONS OF THE NON LOCAL SUPERCURRENT 



The average of the non local supercurrent vanishes because of disorder averaging in the diffusive system*^ or 
because of averaging over the Fermi oscillations in the ballistic system. The disorder average of the square of the non 

local supercurrent (J^ 110 " 100 ^ i nvo lves the diagrams on Fig. 02-(c), the average over disorder of which does not decay 
exponentially over the elastic mean free path Id. 

The non local supercurrent l|13(l can be recast in the form 

it nloc) = A y: £(<T#+*S#$). (21) 

a,fc,a' ,6' cr,r 

where A = Tt(e/h)k\t a ^ a \ 2 \t b ^\ 2 , and X^^, b , = fa^fa^p'g^b ^vy^ ancl wner e er =T, I is the projection of the spin 
along the z axis, and r = 1, 2 is the Nambu index. The notation r in the definition of I^ T j b , corresponds to r = 1 
if t = 2, and f = 2 if r = 1. We deduce from Eq. f2*T|) 

F 1 ~~ ~ (22 ) 



^(nonloc)\ £ ^ ^ Re 



,,b,a',b' ct,t a' 



^a,b,a' ,b'^a,b,a' ,b' "+" ^~a,b,a' ,b'^a,b,a' ,b' 



corresponding to the diagrams on Fig. |3-(c). The disorder averages in each electrode are then carried out, following 
the Appendix. Factoring out the propagators in the superconductor—, and carrying out the summation over the 
conduction channels, leads to 



(/nonloc^ = ^(e/hfA 2 \t a J%^N> h flp f*,^ 

x l exp (- f) + vm^m exp ( ]> sin ^ * (23) 




where we discarded the terms decaying exponentially over the Fermi wave-length. The product fapifLp) 2 ^ s P r0 " 

portional to 1/D' 2 , with D 2 of order N cb a 2 {D /r) 2 , so that (/g onloc ) 2 scales like N ch (see Fig. El for the notations D 
and r). The variance of the supercurrent given by Eq. ll'-'il) involves a "long range" contribution decaying over 1 V1 
and a short range contribution decaying over The former propagates over a much larger distance than the latter. 
Both contributions are identical for normal metals, but the long range contribution dominates for ferromagnets with 

ih<R<i<p. 



V. CONCLUSIONS 



To conclude, we have investigated the possibility of coupling coherently two superconductors by two spatially 
separated ferromagnets. The statistical fluctuations of the Josephson current are proportional to the square root of 
the surface of the junctions. The fluctuating part of the Josephson current is "long range" in the sense that it does 
not decay over the exchange length but decays over the phase coherence length l v . We predict a Josephson current 
mediated by fluctuations if the length of the ferromagnets is larger than ^ and smaller than l v . The supercurrent 
is expected to fluctuate as a function of the relative spin orientation of the ferromagnets. This effect can be used as 
a test of the phase coherence of crossed Andreev reflection, without the competition between the crossed Andreev 
reflection and elastic cotunneling channels 4 ' 7 since the Josephson effect probes solely the anomalous propagator in the 
superconductor. 

Another proposal has been made recently to probe a long range Josephson effect in ferromagnets with non collinear 
magnetizations in a single SFS junction, generating triplet correlations^^ that can also propagate up to l v . This 
effect is not equivalent to the one considered here since it involves propagation in a single electrode. The fluctuations 
of the Josephson current discussed here are also not equivalent to universal conductance fluctuations^ since the root 
mean square of the supercurrent distribution is proportional to the square root of the junction area. 
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APPENDIX A: DETAILS OF DISORDER AVERAGING 
1. Elastic scattering time 

The Green's function G*- -* (k, to) of a disordered isolated ferromagnet, diagonal in the Nambu representation, is given 
by G(°)(k,w) = g(k,u) i + ±M 0) (k,u) , where g{k, lu) is the Green's function of a ballistic isolated ferromagnet, 
and 



(Al) 



is the disorder self-energy, where n is the concentration of impurities. Evaluating the integral in Eq. (|A1|) by contour 
integration leads to 



<3r(°M.M(k, w ) 
G (0),T,2^ (k;W) 



w-h-^k- i/n,i 
1 



with 



r i,i 



T~2.2 



1 + uo/2e F - h/2e F 

T 

1 - uj/2e F + h/2e F '' 



(A2) 
(A3) 

(A4) 
(A5) 



with r = Airc F I (kpv 2 ) the elastic scattering time with h = u> = 0. 



2. Disorder averaging of the supercurrent 



The supercurrent involves the disorder average g^' b ' {ui)g b a ' (w), evaluated in Fourier space in the ladder 
approximation 33 . Using contour integration, we find 



dk 

J2nf< 



G T ' 1 ' A (k,w)G T ' ;i ' A (k + q,w) 



T,2,A/ 



-l + l r(u;-h + l v F /l^) + q2T2e2F 



3k p 



After the summing the ladder diagrams, the poles are found at wave- vector 



7 (diff) 



\ 



ih 



£(ball) 



(A6) 



(A7) 



where we replaced u> by A, as obtained in the energy integration of the local supercurrent [see Eq. JTTJ]. Eq. (|A7|) 
leads directly to Eq. (|19(1 if one assumes the ballistic superconducting coherence length £ ( bal1 ) of the order of 1 /im 
to be large compared to the ballistic exchange length v F /h and to the ferromagnet ballistic phase coherence length 
^(baii)^ rpj^ g eome t r i ca i prefactor (irpp) 2 /kpl e d a ^ is then obtained by evaluating the residue in the integral over q. 



3. Summation over the conduction channels 



The supercurrent given by Eq. involves a summation over the conduction channels. This summation is evaluated 
through 



d~7 exp 



iKh d, 



N, 



cli 



2nydy 



= ^ hT % 7 exp 
1 - iK h ^ 



■. exp 



\ - iK h ) V^ 2 + V 2 



iK h R 



(A8) 



where d a .b = V ' R 2 + V 2 is the distance between the sites a and 6. 
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4. Disorder averaging from real space Green's functions 

The product of the ballistic Green's functions g] l 'l' A (ui)gl'^' A (uj) is proportional to exp [— q^ bal ^d a ^)], with g( ba11 ) = 
iAk + 2iuj/vp + 2/l<p . The wave- vector g( dlff ) given by Eq. (|A7|) is related to g( ba11 ) according to 
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